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Abstract. We discuss the possibility to observe hadron modification in hot 
matter via the correlation of identical particles. We find that the hadron mod- 
ification causes a new type of correlation, back-to-back correlation. 



1. Introduction 

The Hanbury-Brown Twiss (HBT) effect has been widely measured in heavy ion 
collisions. It has been expected that the effect will give some clue to the size 
of the system at freeze-out. Another interesting topic in heavy ion physics is the 
possibility of hadron modification in medium. So far these two have been considered 
as two different aspects of heavy ion physics. The HBT effect is concerned with 
freeze-out where interaction disappears, whereas the hadron modification is caused 
by interaction. However, this picture is purely classical. Since the HBT effect is of 
quantum nature, we need quantum mechanical consideration. 

In relativistic heavy ion collisions, freeze-out looks rather prompt. The distri- 
bution of final state hadrons is, in most of the cases, almost exponential ||^ and this 
suggests that the system is almost thermalized up to some time and then breaks 
up suddenly |2|. Motivated by this, we have modeled freeze-out as follows. The 
system remains thermalized until freeze-out. Hadrons are modified due to interac- 
tion and their masses are shifted. As a result, it is dressed pseudo-particles that are 
thermalized. Then, freeze-out is assumed to take place suddenly. In this paper, we 
investigate how this picture of freeze-out is described quantum mechanically and 
how this scenario is ascertained through the HBT correlation. In Section 2, we 
present the theoretical framework of our work. In Section 3, we show numerical 
results and give some discussion. In Section 4, we conclude and give a perspective. 
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2. Calculation of Correlation Functions 

In this paper, we use the following scalar theory and we ignore the isospin structure, 
because it is not essential in our argument. Throughout this paper, we adopt the 
mean field approximation for simplicity and for clarity. In addition, we do not take 
into account the finiteness of the system. Calculations for finite systems will be 
presented elsewhere 

The theory in the vacuum is given by the following free Lagrangian 

/:o^ld,<j>ix)d^'c^{x)-^mlc^^x), (1) 

where mp is the vacuum mass of the scalar field 4>{x). The Hamiltonian of this 
system is, of course, given by 

i I {n^x) +^^l^'i^))d'x. (2) 

The canonical quantization of this field is normally carried out for the canonical 
variable 0(k) in momentum space and its canonical conjugate 7r(k) defined by 

•/"(k) = j^^^ J(f>{x)exp{~ik-K)d^x, 



(27r) 
1 

(2^3/2 

by introducing an annihilation operator a(k) and a creation operator a^(k) as fol- 
lows: 

0(k) = -^(a(k)+at(-k)), 



r(k) = -z J^(a(k)-at(-k)), (4) 



where o^k = \rm^+^ and we have suppressed the time dependence of the opera- 
tors. After normal ordering, we get the well-known diagonalized Hamiltonian, 

Ho= [ Lui,a\k)a{k)d^k. (5) 



When the temperature and/or chemical density is non- vanishing and the mass of 
the (j) field is shifted, within the mean field approximation, this is expressed by the 
following Lagrangian in medium Cm- 

Cm = ^d^^{x)d^^{x) ~^{ml + ml)c^^x). (6) 
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The mass shift SM is given by 

SM ^ ^ml + ml- 7710. (7) 

Let us define the quanta that diagonalize the Lagrangian (^) as 6(k). The point 
here is that the fe-quanta are, in general, different from the a-quanta which are the 
fundamental excitations in the vacuum. In other words, the a-operators do not 
diagonalize the Hamiltonian in medium Hm- I3y inserting Eq. (^) into Hm and 
carrying out normal ordering, we obtain 

Hm = Ho + Hi, 

Hi ^ / — [a(k)a(-k)+2at(-k)a(-k)+at(k)at(_k)]d3fc. (8) 

Therefore, in medium, mode k and mode — k of the a-quanta are mixed. This 
Hamiltonian (|^) can be exactly diagonalized with the following Bogoliubov trans- 
formation: 

a^(k) = coshrkfe^(k) + sinhrk k), a^(— k) = coshr_k fi^l— k) + sinhr_k f'(k), 
a(k) = sinhrk k) -I- coshrk fc(k), a(— k) = sinhr_k ^^(k) + coshr_k k), 

(9) 

where rk is a real number which satisfies the following equation, 

2 

tanh2rk = (10) 
2u;^ -I- mf 

This gives the exact relationship between the quanta in the vacuum (a-quanta) 
and the quanta in medium (6-quanta) in this theory. With the 6-operators, the 
in-medium Hamiltonian (^) is diagonalized as 

Ho = j ni,b'<{k)b{k)d^k, (11) 

where fik is given by 



= Jml + ml+k.-^. (12) 



This is almost trivial, since the mass of the (p field is shifted to -y/ rriQ + m\ in 
medium. However, the important point here is that it is the 6-quanta that are 
thermalized in medium, since the ^-operators diagonalize the in-medium Hamilto- 
nian. On the other hand, the ^-quanta are not observed experimentally. It is the 
a-quanta that are observed. Therefore, in calculating final state observables, we 
have to evaluate the expectation value of operators defined in terms of a and 
operators, ©(a, a^) with the density matrix defined in the 6-basis, ph, i.e., 

{0{a,a^)) ^TipbO{a,a)). (13) 
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The calculation of the one and two-particle distribution functions which are needed 
to obtain the two-particle correlation is straightforward, but the Glauber - Sudar- 
shan representation of the thermal density matrix , 



Pb 



1. *^ 



is useful 1^. Here |/3(k)) is a coherent state satisfying fe(k)|/3(k)) = /3(k)|/3(k)) and 
P(/3(k)) is defined by 



where Jb (k) is the Bose-Einstein distribution function with mass a/ ttIq -I- rnf . With 
this formula, we obtain the one-particle distribution in the final state, 

{a\k)a(k)) = cosh2rk /^(k) + sinhVk- (16) 

Since we have ignored the finite size effect, the two-particle distribution function 
(at(k)at(k')a(k)a(k')) takes a trivial value (a^(k)a(k))(at(k')a(k')) unless k = 
±k'. For k = ±k' cases, we obtain 

(a^(k)a1'(k)a(k)a(k)) = 2(at(k)a(k))2, 
(at(k)at(_k)a(k)a(-k)) = (at(k)a(k))2 + /s(k)(/B(k) + l) sinh22rk + 2sinhVk, 

(17) 

where we have assumed the uniformity of the system, i.e., /^(k) = /b(— k) and 
Tk = "T-k- From this, we get the two-particle correlation fmictions C2(k, ±k) as 
follows: 

C2(k,k) = 2, 

^ „ , , , , /B(k)(/B(k) + 1) sinh^2rk + 2sinhVk 

G2(k, — k) = IH 2 

(cosh2rk /s_B(k) -I- sinh rk)^ 

^ 1. (18) 

In all other cases, the two-particle correlation function is 1. This result is quite 
interesting. It implies two important issues, (i) the intercept of the correlation 
function C2 (k, k) remains the canonical value even if quanta in medium are different 
from those in the vacuum (ii) back-to-back correlation (C2(k, — k) 7^ 1) is generated 
by hadron modification. This is caused by the mixing of k and — k modes due to 
the mean field effect. In our case, the mean field carries no momentum. Due to 
momentum conservation, no other mixing is possible. Hence, only k and — k modes 
get additional correlation. 
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3. Results and Discussion 



The medium effects on the momentum distribution of pions and kaons are evaluated 
numerically for static fireballs in Figs. 1 and 2. These figures indicate that the 
medium modification of hadronic masses results in a very small change in the shape 
of the momentum distributions. Fig. 3 indicates the momentum dependence of 
C2(k, — k) for pions, assuming a freeze-out temperature of Tf = 120 MeV and less 
than 1 % values for the mass-shift of pions in medium. Figure 3 indicates that 
the deviation of C2(k, — k) from unity is very small for pions and that the effect is 
smaller than the current levels of error bars for the two-particle correlation function. 
However, the medium modification of kaons may be observable in the two-particle 
correlation function according to Figure 4. The enhancement at large values of 
particle momenta in the CMS of the fireball is a strongly increasing function of the 
medium modification of the kaonic masses. The effect is decreasing with increasing 
values of the freeze-out temperature Tf. 

Our result is essentially different from that of Ref. |^ where over-bunching of 
identical pions was predicted at small relative momenta, i.e. C2(k, k) > 2, while 
in our case the intercept parameter takes up the canonical thermal value of 2. We 
find, in contrast to Ref. that mode- mode mixing generates positive correlation 
of kaons at large relative momenta, if the mass of the kaons is modified in the 
medium. 




Fig. 1. Solid line indicates the spectrum of pions from a fireball, when a small 
in-medium mass-shift of -10 MeV is invoked, dotted line stands for the no 
mass-shift case. 
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Fig. 2. Solid line indicates the spectrum of kaons from a fireball, when a 
small in-medium mass-shift of -14 MeV is invoked, dotted line stands for the 
no mass-shift case. 
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Fig. 3. Solid line indicates the correlation function C2(k, — k) of pions from a 
fireball, when a small in-medium mass-shift of -10 MeV is invoked, dashed line 
stands for the 5M = —5 MeV, dotted line stands for the no mass-shift case. 
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Fig. 4. Illustration of the dependence of C2(k, — k) on the value of the 
freeze-out temperature 7/, for kaons. 
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Fig. 5. Same as Fig. 3, but for kaons, and a little larger mass-shifts. 
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4. Summary 

A new method has been found to test the medium modification of kaons, utihzing 
their quantum correlations at large momentum difi^erence. The effect follows from 
basic principles of statistical physics and canonical quantization. 

Back-to-back correlations are not contaminated by resonance decays, in con- 
trast to the Bose-Einstein correlation functions at small relative momenta. For a 
locally thermalized expanding source, however, back- to-back correlations will ap- 
pear in the rest frame of each fluid clement. Thus, for such systems a more realistic 
estimate is necessary to learn more about the magnitude of the correlations of kaons 
at large relative momenta. 

In the present study we have neglected finite size effects, which will make the 
correlation function vary smoothly around both ki = k2 and ki = — k2. 
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